Measurement of the quasi-elastic axial vector mass in neutrino-oxygen interactions 
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The weak nucleon axial- vector form factor for quasi-elastic interactions is determined using neu- 
trino interaction data from the K2K Scintillating Fiber detector in the neutrino beam at KEK. 
More than 12,000 events are analyzed, of which half are charged-current quasi-elastic interactions 
v^n — * fi~p occurring primarily in oxygen nuclei. We use a relativistic Fermi gas model for oxygen 
and assume the form factor is approximately a dipole with one parameter, the axial vector mass 
Ma, and fit to the shape of the distribution of the square of the momentum transfer from the 
nucleon to the nucleus. Our best fit result for Ma = 1-20 ± 0.12 GeV. Furthermore, this analysis 
includes updated vector form factors from recent electron scattering experiments and a discussion 
of the effects of the nucleon momentum on the shape of the fitted distributions. 

PACS numbers: 13.15.+g; 23.40.Bw; 25.30.Pt 



INTRODUCTION 



II. CROSS SECTION AND FORM FACTOR 
EXPRESSIONS 



The structure of the nucleon, as measured both by elec- 
trons and neutrinos, has been a subject of experimental 
study for decades. The discovery of neutrino oscillation 
and the availability of high precision electron scattering 
measurements have renewed interest in the study of neu- 
trino interactions on nuclei. Neutrinos offer unique in- 
formation about the structure of the nucleon and the nu- 
cleus. There are many experimental neutrino programs 
now running, under construction, or being planned for 
the near future, all of which use nuclear targets such 
as oxygen, carbon, aluminum, argon, or iron. Likewise, 
there has been significant progress in the calculation of 
cross sections, nuclear corrections, and backgrounds to 
specific processes. Improvement of these models, sup- 
ported by neutrino data, will be important for the up- 
coming precision neutrino oscillation studies. 

In this study we analyze distributions of the square of 
the four-momentum transfer Q 2 = —q 2 = — (p^ —p v ) 2 re- 
constructed from neutrino-oxygen interactions, where p M 
and p„ are the momenta for the outgoing muon and inci- 
dent neutrino. Using data from the Scintillating Fiber 
(SciFi) detector in the KEK accelerator to Kamioka 
(K2K) neutrino beam, we fit for the value of the axial 
vector mass Ma, the single parameter in the axial vec- 
tor form factor (assuming a dipole form) for quasi-elastic 
(QE) interactions. For QE interactions, this parameter 
is obtained only from neutrino-nucleus scattering exper- 
iments. This is the first such measurement for oxygen 
nuclei, and we include a discussion of the effects of the 
oxygen nucleus and nucleon momentum distribution on 
the shape of the Q 2 distribution. 

In the next section we briefly discuss the quasi-elastic 
form factors as well as list the cross-sections for non- 
quasi-elastic processes and nuclear effects. Then a section 
describes the K2K experiment, the neutrino beam, and 
the SciFi detector, including the detector performance. 
Following that are sections on the analysis technique, and 
the results, which include detailed discussion of the major 
systematic effects. 



A. Quasi-elastic cross section 

The differential cross section da/dq 2 for neutrino 
quasi-elastic scattering (v^ n — > fj," p) is described in 
terms of the vector, axial- vector, and pseudo-scalar form 
factors. The differential cross sectionjlj is written as: 



da" 



M 2 G 2 F cos 2 1 
8ttE 2 
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where, s and u are Mandelstam variables, (s-u) = 
AME V + q 2 — m 2 , m is the outgoing lepton mass, M is 
the target nucleon mass, and E v is the neutrino energy. 
A(q 2 ), B(q 2 ), and C(q 2 ) are: 
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*Now at University of Minnesota, Duluth 



In these expressions, the pseudo-scalar form factor Fp 
is negligible for muon neutrino scattering away from the 
muon production threshold and is not included. Fa is 
the axial vector form factor we will extract from the data. 
Fy(q 2 ) and F v (q 2 ) are the Dirac electromagnetic isovec- 
tor form factor and the Pauli electromagnetic isovector 
form factor, respectively. These formulas also assume the 
conserved vector current (CVC) hypothesis, which allows 
us to write F v and F v in terms of the well measured 
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Sachs form factors G E , G E , and G^: 



fUv 2 ) = (i- 
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In this paper we use the updated measurements of the 
Sachs form factors from 0, • These new form factors 
have a significant effect on the extraction of Fa, com- 
pared to the previous dipole approximations. For the 
range of Q 2 of interest in this experiment, the updated 
values differ from the old form factors by up to ± 10%. 
We present results with both the new and the old form 
factors in this paper. 

We approximate the axial vector form factor Fa as a 
dipole 



F A (q 2 



1.2720 



(1 - (qt/MDY 



(4) 



which has a single free parameter, the axial vector 
mass Ma- Previous studies show that this approxima- 
tion is reasonable [4J, |5j, li| The constant FA(q 2 —0) = 
9a/ 9v = 1.2720±0.0018 is determined from neutron de- 
cay measurements p|. Because the Sachs form factors and 
other constants are precisely measured, the single param- 
eter Ma can be determined from quasi-elastic neutrino 
interaction data. 



B. Other cross sections 

For this analysis, approximately half of the data comes 
from non quasi-elastic interactions, especially single pion 
events from the production and decay of the N* and 
A baryon resonances within the nucleus. This back- 
ground is described by the NEUT neutrino interaction 
Monte Carlo simulation || used by the K2K and Super- 
Kamiokande experiments. The resonance single pion 
events are from the model of Rein and Sehgal|9j. 

Additional backgrounds, less important for our beam 
energy around 1.2 GeV, are also included. Deep inelas- 
tic scattering is from GRV94[lO| for the nuclear struc- 
ture functions with a correction described by Bodek and 
Yangjll]]. The software PYTHIA/JetSet 12] is used to 
generate these events. This analysis takes the charged 
current coherent pion cross section to be zero following 
[l3j . and include neutral current coherent pion interac- 
tions as in Rein and Sehgal 14] . 



C. Nuclear Effects 

Equation 2] is the differential cross section for the free 
nucleon, and must be modified to account for the effects 



of a nucleon bound in a nucleus. In the SciFi detector, 
the fiducial mass fractions are 0.700 H 2 0, 0.218 Al, 0.082 
HC, with an error of ±0.004. Our neutrino interaction 
Monte Carlo treats the entire fiducial mass as if it was 
made of H2O; for targets other than a proton in hydro- 
gen, we use a uniform Fermi gas model with kf — 225 
MeV/c for the nucleon momentum and an effective bind- 
ing energy of -27 MeV, which is appropriate for oxygen. 
The primary effect of this nucleon momentum distribu- 
tion on the quasi-elastic events is a significant suppres- 
sion at low Q 2 due to Pauli blocking and a smaller overall 
suppression of ~2% for the entire Q 2 distribution. The 
Fermi gas model is also applied to the non quasi-elastic 
interactions. 

In addition to cross section effects, there are final-state 
interactions. The nucleus will cause reinteraction or ab- 
sorption of secondary pions and recoil protons and neu- 
trons in the neutrino interaction final state. Our model 
for these reinteractions is described in [8j with references. 
This will affect the observed distribution of the number of 
tracks. The resulting fi~ is also affected by the Coulomb 
interaction as it leaves the nucleus, losing approximately 
3 MeV, though this effect is implicitly included in the 
Fermi gas parameters. The above nuclear effects are dis- 
cussed quantitatively in the results section. 



III. EXPERIMENT 
A. The beam and detectors 

The KEK to Kamioka (K2K)[H E E3 experiment 
is a long baseline neutrino oscillation measurement in 
which a beam of neutrinos is sent from the KEK ac- 
celerator in Tsukuba, Japan to the underground Super- 
Kamiokande detector[18|. The neutrinos pass through a 
set of near neutrino detectors 300 meters from the target, 
after which they travel 250 km to Super-Kamiokande. 
The analysis in this paper considers only neutrino inter- 
actions detected in the Scintillating Fiber (SciFi) detec- 
tor, one of the near detectors. 

The wide-band neutrino beam at KEK is produced 
when 12 GeV protons hit an aluminum target. Two mag- 
netic horns focus positively charged pions and kaons into 
a 200 meter long decay pipe, where they decay to fi + and 
Vp. The are absorbed by the beam dump plus ap- 
proximately 100 meters of earth between the decay pipe 
and near detector hall. The resulting neutrino energy 
is between 0.3 and 5 GeV and peaks at 1.2 GeV. The 
contamination in this beam includes 1.3% v e and 0.5% 
anti-i^,, estimated from a Monte Carlo simulation of the 
beam. 

The near detector hall of the K2K experiment contains 
several detectors, shown in Fig. ^ The first detector 
in the beam is the one-kiloton water Cerenkov detector. 
This study uses data from the SciFi detector, which is 
described in detail below. Following SciFi is the location 
of the lead glass detector which was used to measure the 



4 



v e contamination in the beam. The lead glass detec- 
tor was removed in 2002 and in its place was a proto- 
type for a plastic scintillator (SciBar) detector. Then in 
2003, the full SciBar detector was installed, though 
data from this last running period is not used in the 
present analysis. Finally, there is a muon range detector 
(MRD) |2(j which is used to estimate the momentum of 
the muons which escape the SciFi detector from charged 
current neutrino interactions. The MRD is also used to 
monitor the stability of the neutrino beam. 



Veto/Trigger 
Muon Range Detector counter 




FIG. 1: The arrangement of the near neutrino detectors at 
KEK. The beam comes in from the right and continues to 
Super-Kamiokande, 250 km away to the left. The SciFi de- 
tector in on a stand in the middle. 



The prediction for the shape of the neutrino energy 
spectrum of the K2K beam has significant uncertainty, 
up to 20% at higher energies. This prediction uses 
a Sanford-Wang parameterization of hadron production 
data and is verified using pion monitors downstream from 
the target p^. This energy spectrum is measured using 
data from the near detectors and is used as input to the 
oscillation analysis 0|- The energy spectrum analysis 
is a simultaneous fit to the muon momentum and muon 
angle distributions from charged current interactions in 
the one-kiloton water Cerenkov detector, the SciFi de- 
tector, and the SciBar detector. The free parameters 
in this fit are a scale factor for the flux in eight energy 
regions, a scale factor for non quasi-elastic events, and 
many systematic error parameters specific to each detec- 
tor. For this paper, we will refer to the above procedure 
as the neutrino energy spectrum measurement, and it 
defines the baseline Monte Carlo prediction for the SciFi 
detector data, prior to any fitting for the axial-vector 
mass, and is used throughout the discussion and plots in 
this Sec. III. This default MC simulation also uses zero 
charged current coherent pion andM^ = 1.1 GeV. The 
resonance single-pion cross section also involves its own 
axial- vector term with its own — 1.1 GeV. The anal- 
ysis described in Sec. IV and V is mostly independent 
from the energy spectrum analysis, but uses a similar 



strategy. 

B. The SciFi detector 

The SciFi detector [21II22] ] consists of scintillating fiber 
tracking layers between aluminum tanks filled with wa- 
ter. A schematic diagram is included in Fig. [3 There 
are a total of twenty 240 cm x 240 cm wide tracking lay- 
ers, each of which consists of fibers oriented to give the 
particle location in the horizontal and vertical direction. 
These fibers are glued, one layer on each side, to a hon- 
eycomb panel which is 260 cm square. The distance be- 
tween two tracking layers is 9 cm. Between the first and 
the twentieth layer are nineteen layers of aluminum tanks 
whose walls are 0.18 cm thick with an interior thickness 
6 cm filled with water. 




FIG. 2: A schematic diagram of the SciFi detector. 



The scintillating fibers have a diameter of 0.7 mm and 
are read out by coupling them to image intensifier tubes 
and CCD cameras. The image intensifier preserves the 
position information of the original photo-electron. At 
the final stage, the light is recorded by a CCD camera. 
A total of 24 of these are used to read out 274,080 scin- 
tillating fibers. To reconstruct which fibers were hit, a 
one-to-one correspondence between the fibers and the po- 
sition of pixels on the CCD camera is obtained from pe- 
riodic calibration using an electro-luminescent plate. 

To select charged-current events for this analysis, we 
require at least one track start in the SciFi fiducial vol- 
ume and extend to, and stop within, the MRD. The 
fiducial volume has a mass of 5590 kg, includes the first 
through 17th tanks of water, and the reconstructed ver- 
tex must be within 110 cm from the center of SciFi in 
horizontal and vertical directions. 

This requirement means that all events selected for 
this analysis have hits in at least three tracking layers of 
SciFi. There are also upstream and downstream scintil- 
lator hodoscopes which are read out by photo-multiplier 
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tubes; we require a matching hit downstream and no hit 
upstream. This allows us to veto the small number of 
muons which come from in-time muon generation in the 
rock and upstream material in the detector hall, and also 
muons surviving from the target. Cosmic ray muons are 
suppressed by the beam timing requirement and are also 
negligible. 

Tracks are reconstructed in the horizontal and vertical 
projections separately and then matched. The efficiency 
for reconstructing muon tracks with hits in three SciFi 
tracking layers is ~70%, and rises to nearly 100% for 
tracks that penetrate five or more layers. Second tracks 
are required to produce hits in at least three SciFi layers, 
but there is no restriction on the maximum length. 

Prior to 2002, all muons from SciFi are required to pass 
through and produce hits in segments of the lead glass 
detector, and these segments must match the location of 
the reconstructed track seen in SciFi and the MRD. On 
average, they deposit around 0.4 GeV of energy in the 
lead glass, though only path length, and not pulse size, 
is used to estimate the energy loss in this case. In the 
K2K-II run period, muons traveling through the SciBar 
prototype lose around 0.023 GeV of energy, though we 
do not require the track to pass through this detector. 
Muons traveling through many layers in SciFi deposit up 
to 0.3 GeV of energy. 

The Muon Range Detector is made of alternating lay- 
ers of drift tubes and iron plates; the first detection layer 
is upstream of the first piece of iron. The first four layers 
have a thickness equivalent to about 0.14 GeV of energy 
loss each, and the remaining layers are twice as thick. 
The muon momentum can then be estimated by calcu- 
lating the muon's range from the interaction vertex. 



C. Data samples 

The data for this analysis are obtained from two run- 
ning periods between November 1999 and June 2003. 
The primary distinction between them is the configu- 
ration of the Super-Kamiokande detector, though there 
were simultaneous changes in the near detector configu- 
ration important for this analysis. We refer to the first as 
the "K2K-I" period; muons from neutrino interactions in 
SciFi pass through the lead glass detector on their way to 
the MRD. For these data, we accept muons which pene- 
trate as little as one MRD detection layer, which corre- 
sponds to a muon momentum threshold of 675 MeV/c. 
The second running period is called "K2K-IIa" and has 
the prototype for the plastic scintillator detector SciBar 
[T^| in place of the lead glass. For K2K-IIa, we require 
that the muons produce hits in the first two layers of the 
MRD, which gives a threshold of 550 MeV/c, in order to 
reduce the contamination from pions reaching the MRD. 
Data from the continuation of the K2K-II period are not 
used in this analysis. In all cases, we require the muon 
not exit the MRD, which results in a maximum muon 
momentum of 3.5 GeV/c. 



When two tracks reach the MRD, the longest, most 
penetrating track is assumed to be the muon. We have 
estimated using the MC that approximately 2% of these 
longest tracks are not the muon track, and another 0.5% 
of events were from neutral current interactions which 
had no muon at all. 

This analysis uses only one-track and two-track events. 
Since quasi-elastic interactions will not produce such 
events, the 3% of events with three or more reconstructed 
tracks are discarded. For one-track events, the recoil pro- 
ton or a pion is absent or below threshold. The require- 
ment of three layers for the second track corresponds to 
a threshold of 600 MeV/c proton momentum and 200 
MeV/c pion momentum. 

The MC simulation includes the rescattering of pro- 
tons, neutrons, pions, and other hadrons from the neu- 
trino interaction final state as they leave the nucleus. The 
models for these final state interactions give good agree- 
ment with the number of tracks seen in SciFi, shown in 
Tab[I] Our estimate of the uncertainty in the number of 
two-track events due to the efficiency for finding the sec- 
ond track is 5%, which translates into ± 100 events for 
the K2K-I two-track sample. These events primarily mi- 
grate to or from the one-track sample. The extreme case 
of zero nuclear final state interactions in the neutrino in- 
teraction MC leads to 20% more events in the two-track 
category [23j. Further discussion of these effects can be 
found in Sec. V of this paper. 



sample 


1-track 2-track 3-track 


K2K-I Data 


5933 


2181 


187 


K2K-I MC 


5920 


2176 


203 


K2K-IIa Data 


3651 


1344 


148 


K2K-IIa MC 


3583 


1396 


136 



TABLE I: Comparison between data and MC simulation of 
the number of reconstructed tracks observed in the SciFi de- 
tector for the K2K-I and K2K-IIa data samples. The MC is 
normalized to have the same number of events as the data. 
The estimated error in the number of two-track events due 
to tracking efficiency in reconstructing short, second tracks is 
5%. 

For two-track events, we separate quasi-elastic from 
non quasi-elastic events. Quasi-elastic interactions are a 
two-particle scattering process; the measurement of the 
muon momentum and angle is sufficient information to 
predict the angle of the recoil proton. If the measured 
second track agrees with this prediction within 25°, it is 
likely a QE event. If it disagrees, then it becomes a part 
of the non-QE sample. This is shown in Fig.0 where the 
inset diagram demonstrates the kinematic quantity A9 = 
angle between predicted and measured second track angle 
with respect to the beam. The quantity A8 is plotted in 
this figure with the data and the baseline MC normalized 
to the data. 

The value for this A8 cut is chosen to give good sep- 
aration between the QE and nonQE enhanced samples. 
We used the Monte Carlo simulation to estimate the ef- 
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A6 (degrees) 



FIG. 3: The distribution of A8, the difference between the 
predicted second track angle and the observed angle for K2K-I 
data. The histogram shows the Monte Carlo prediction, while 
the hatched region shows the QE fraction. In this analysis, 
A9 = 25 degrees is used to separate a QE enhanced sample. 
The inset diagram shows the definition of A9. 



ficiency for detecting the QE events, after all the cuts 
described above. Also we have estimated the purity of 
each sub-sample. These are shown in Tab.lTTl After these 
cuts, the total number of events in each sample is given 

in Tab. nm 





1-track 


2-track 


Total 






QE nonQE 




K2K-I 


35 (63) 


5 (63) 2 (17) 


42 


K2K-IIa 


38 (61) 


5 (61) 2 (15) 


45 



TABLE II: Total reconstruction efficiency [%] for quasi-elastic 
interactions in each data set, the portion of efficiency from 
each sub-sample, and the QE purity of each sample (in paren- 
thesis, [%]), estimated with the MC simulation. 





K2K-I 
Q 2 > 0.0 Q 2 > 0.2 


K2K-IIa 
Q 2 > 0.0 Q 2 > 0.2 


1 track 


5933 


2864 


3623 


1659 


2 track QE 


740 


657 


451 


388 


2 track nonQE 


1441 


789 


893 


478 


Total 


8114 


4310 


4967 


2525 



TABLE III: Number of events in three event samples and 
two data periods for the SciFi detector. Only events with 
reconstructed Q 2 > 0.2 (GeV/c) 2 are used for this Ma mea- 
surement, and are shown in separate columns and described 
in Sec. IV. 



D. Muon momentum and angle distributions 

An example of the muon momentum and muon angle 
distributions for the K2K-I data along with the Monte 



Carlo prediction are shown in Fig. g] and Fig. [5] The MC 
distribution is normalized to the same number of events. 




0.5 1 1.5 2 2.5 3 3.5 

muon momentum (GeV/c) 

FIG. 4: Muon momentum distribution for all K2K-1 one-track 
and two-track events. The QE fraction, estimated from the 
MC simulation, is shown as the shaded region. The errors on 
the data are statistical only. 



We observed a deficit of events whose muon is at angles 
near the direction of the beam compared to our Monte 
Carlo simulation; this is also discussed in [13 . The dis- 

700 1= 

600- Jf^ftf, 

500 - j_T \, 

400 - J 4= 

300- ^J - ' -1 ^ L ~ Ll_Ll_ ^^^^^^^^^ 

°0 10 20 30 40 50 60 

muon angle (degrees) 

FIG. 5: Muon angle distribution for all K2K-I one-track and 
two-track events. The QE fraction, estimated from the MC 
simulation, is shown as the shaded region. Only statistical 
errors are shown. 

crepancy was observed in all K2K near detectors, includ- 
ing SciFi, and is presumed to be from some aspect of the 
neutrino interaction model. The analysis of data from 
the SciBar detector |l3| indicated that most, if not all 
of this deficit is because there is too much CC coherent 
pion production in the MC. The SciBar data are consis- 
tent with zero charged-current coherent pion. Examples 
of the disagreement from SciFi data are shown in Fig. [U 
with and without charged current coherent pion events. 
The MC distribution with zero CC coherent pion is nor- 
malized to the same number of events. In this analysis, 
we assume there is zero CC coherent n production. 
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FIG. 6: Distribution of cos($^) for data and MC, showing the data at small angle, which is also the low Q 2 . Left to right are 
the one-track, two-track QE enhanced, and two-track nonQE enhanced samples. The top line is with charged-current coherent 
pion, the second line without. The shaded region is the QE fraction, estimated from the MC simulation. The distributions are 
normalized to the total number of events when all three samples are combined. 



IV. ANALYSIS 

A. Calculating Q 2 and E v 

The kinematics of the muon candidate, the longest 
track in our events, are sufficient to estimate the energy 
of the neutrino E r v ec and the square of the momentum 
transfer Q 2 ec , if the interaction is quasi-elastic. 



El 



(mjy + e B )E f2 - (2m N e B + e| + m 2 l )/2 
m N + e B - E^+p^ cos M 



Qlec = -Q 2 = -lEuiEn - Pn cos 6^) + ml 



,(5) 



(0) 



Here, E^ and p^ are the energy and momentum of the 
muon determined from the range, 9^ is the angle relative 
to the beam direction, determined from the hits in the 
SciFi detector. Note that E v appears in the expression 
for Q 2 ec , and here we use E r v ec . The quantity e B = -27 
MeV for oxygen is the effective binding energy parame- 
ter from the Fermi gas model. The masses and 
are for the nucleon and the muon. The resolution for 
Efj, is 0.12 GeV, due mainly to the MRD segmentation. 
The resolution for M is about 1 degree, but there is a 
tail to this distribution. The resulting value for E v res- 
olution (for QE events) is 0.16 GeV and the resolution 
for Q 2 is 0.05 (GeV/c) 2 also with a tail coming from the 
measured angle. Finally, this formula assumes that the 
target neutron inside the nucleus is at rest, ignoring the 
nucleon momentum distribution for the event reconstruc- 
tion. Fluctuations due to Fermi motion are about half the 
size of those due to detector and reconstruction effects, 
and contribute only a small amount to the reconstructed 
energy resolution. 

It is important to note that these formulas are used for 
all events even though half the interactions are not quasi- 
elastic, because we do not identify the interaction mode 
on an event-by-event basis, nor is our beam at a fixed 
energy. The reconstructed E v and Q 2 are systematically 
off for these non quasi-elastic events: E r v ea is low by ^0.4 
GeV and Q 2 rec is low by ~0.05 (GeV/c) 2 . However, all 



events are treated the same way, both data and Monte 
Carlo events. Thus, the comparison of data and MC in 
the fit is valid, but the distributions of the reconstructed 
values are affected by the non quasi-elastic fraction. 



B. Fit procedure 

After calculating ££ ec and Q 2 ec for each event, the 
data are binned in five E r v ec bins: 0.5 to 1.0, 1.0 to 1.5, 
1.5 to 2.0, 2.0 to 2.5, and greater than 2.5 GeV. The 
data are divided into Q 2 bins each of width 0.1 (GeV/c) 2 . 
To ensure there are at least five events in each bin, the 
smaller number of events at higher Q 2 are combined into 
a single bin. 

The expectation for the number of the events in each 
bin is computed from the Monte Carlo simulation for 
different values of the axial-vector mass and some sys- 
tematic error parameters. We perform a maximum like- 
lihood fit to the data by minimizing the negative of the 
logarithm of the likelihood which is based on Poisson 
statistics for each bin. In our case we use the modified 
form given in the Review of Particle Physics 

N 

- 21nA(0) - 2 Y^WM -rii+rii lnfa/i/^fl))] (7) 
»=i 

in which fi(6) and n, are the predicted and observed 
values in the z-th bin for some values of the parameters 
9. The minimum of this function follows a chi-square 
distribution and can be used to estimate the goodness of 
the fit. 

The expectation for each reconstructed E v and Q 2 bin 
is computed as follows: 

N t otai(ntra,ck, E r v ec , Q 2 ec ) — A NQE(n tra ck, E^ ec , Q 2 ec ) 



B x N non QE(n tra ck, El, ec ,Q 



2 ) 

recJ 



(8) 



where Nqe and N non QE are the separate contributions 
of quasi-elastic and non quasi-elastic events. 
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The free parameter A is the overall normalization. Five 
parameters $(i?*'' ue ), not included in the above expres- 
sion, are used rescale the neutrino flux in each energy 
region, four of which are unconstrained in the fit, while 
the relative flux for energies from 1.0 GeV to 1.5 GeV 
is fixed at 1.0. The flux is reweighted based on the true 
energy of the MC events, and applies to both QE and 
nonQE events. The nonQE background is reweighted us- 
ing the unconstrained parameter B, which is referred in 
the rest of this paper as the ratio nonQE/QE: the relative 
reweighting of our default MC calculation. Because of the 
separation of the two-track QE and non-QE samples, the 
nonQE / QE ratio will be constrained by the background 
and allow a fit for the QE axial-vector form factor. The 
parameters A and B are relative to the data/MC nor- 
malization calculated using all other parameters at their 
nominal values, including Ma = 1-1 GeV. Importantly, 
changing the value for M9 E changes the absolute cross- 
section for QE, which will in turn affect the fit value for 
the free nonQE/QE parameter and the overall normal- 
ization. 

In this expression, Nqe is based on a calculation of the 
quasi-elastic cross section with the free parameter Ma- 
This cross section is computed using the true energy and 
Q 2 and convoluted with the detailed shape of neutrino 
energy spectrum, flux(E), from the beam MC calcula- 
tions and the hadron production parameterization used 
in [13. 



ill bins 



all bins 

NQE{ntrack-> E r u ec ,Q rec ) = 

Etrue,Q 2 true 

xda/dQ 2 (Et rue ,Q 2 true ,M A ) x i?(£*™ e , Q 2 true ) 

x M(£*« e , Q 2 true -» n track , £- c , Q 2 rec )] . (9) 

Nuclear effects that modify the cross section, especially 
Pauli blocking and and other effects of the nucleon mo- 
mentum distribution, are included using the factor R, 
and are discussed in Sec. V. 

Because the cross section is calculated using true kine- 
matics, it must be modified to account for detector ac- 
ceptance and resolution, as well as nuclear final state 
interactions, in order to obtain the expectation in differ- 
ent reconstructed E v and Q 2 bins. This is done with a 
migration matrix M in the above equation where n tra ck 
refers to the one-track, two-track QE, and two-track non- 
QE samples. This matrix is computed directly from the 
Monte Carlo simulation. This result is then applied to 
the calculated cross section to determine the number of 
QE events in each reconstructed E v and Q 2 bin. In con- 
trast, the shape of the non-QE background is taken di- 
rectly from the Monte Carlo simulation and already in- 
cludes these effects. 

The combination of four flux reweighting factors $(-E„) 
and the overall normalization are unconstrained. The 
parameter Ma itself affects the total cross-section as a 
function of energy. In this way, we are fitting the shape 
of the Q 2 ec distribution separately in each energy region. 



This ensures that the axial mass measurement is not sig- 
nificantly biased by the normalization in any one energy 
bin. 

Finally, the lowest Q 2 ec bins, events below 0.2 
(GeV/c) 2 , are not included in the fit. The low Q 2 re- 
gion is where there is the largest uncertainty due to the 
model for nuclear effects, especially Pauli blocking. This 
eliminates almost half the data, and the total number of 
events actually included in the fit is shown in the second 
column for each data set in Tab. 11111 Low Q 2 ec events 
are also events at low angle, shown by the cos(#) term in 
parenthesis in Eq. and corresponds to the right-most 
two bins in the cos(#) histograms in Fig. for neutrino 
energies around 1.0 GeV. 



C. Fit Parameters 

We fit a large collection of E T v ec and Q 2 ec distributions: 
two data sets K2K-I and K2K-IIa, each with one-track, 
two-track QE, and two-track non-QE subsamples, a to- 
tal of 242 bins. The Monte Carlo predictions for these 
data sets are computed separately using MC samples that 
are more than 15 times larger than the data. The free 
parameters $ for the flux at each energy are common 
to both data sets, as is the overall normalization factor, 
the non-QE/QE ratio and proton rescattering. There are 
separate 2-track to 1-track migration parameters for each 
data set, ten parameters in total. 

These last three parameters are constrained by adding 
a term to the chisquare. A reweighting or migration is 
computed and a systematic error for each of these, as- 
sumed to be approximately Gaussian, is estimated from 
studies of the detector and interaction Monte Carlo simu- 
lations. As the fit is performed, a chisquare is evaluated, 
assuming this Gaussian shape and error, and is added to 
the total chisquare. The total degrees of freedom is thus 
(242 + 3) chisquare terms - 10 parameters = 235 degrees 
of freedom. 

Proton rescattering is taken to be uncertain by ± 20% 
from the nominal value used in the NEUT Monte Carlo 
simulation. A reweighting of events is calculated from a 
second full-detector Monte Carlo simulation using 80% of 
the proton reinteraction cross sections, and the system- 
atic error parameter is used to interpolate between these 
data samples. This reweighting has the primary effect of 
increasing or decreasing the number of QE events in the 
two-track sample, by up to +4.5% for the 80% case. The 
actual effect on the fit is to increase or decrease the num- 
ber of events in the QE enhanced sample because only 
that part of the two-track sample has high QE purity. 

The 2-track to 1-track migration is assigned a ± 5% 
error. This error is based on the estimated error in the 
track-finding efficiency for short, second tracks. Because 
the detector response and tracking is tuned separately for 
the K2K-I and K2K-IIa data samples, we include sepa- 
rate parameters for each sample in our fit. This migration 
is applied to all events in the MC, not just the QE events 
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K2K-IIa one-track Q 1 (OeV/c) ! K2K-IIa two-track QE Q ! (GeV/c) ! K2K-I[a two-track nonQE Q ! (GeV/cf 

FIG. 7: The data and the best fit Q^ec distributions for K2K-1 data (top) and K2K-IIa data (bottom) for the 1-track, 2-track 
QE enhanced, and 2-track non-QE enhanced samples. The shaded region shows the QE fraction of each sample, estimated 
from the MC. The contribution from each energy region is summed for each plot. The lowest two data points in each plot are 
not included in the fit, due to the large uncertainty in the effects of the nucleus. 



as for proton rescattering. This parameter is also corre- 
lated with neutrino interaction effects such as proton and 
pion rescattering. 



V. RESULTS AND DISCUSSION 
A. Fit Results 

The result of the combined fit is M A = 1.20 ± 0.12 
GeV. The chisquare value for this fit is 261 for 235 degrees 
of freedom. The Q 2 distributions for the data and the 
MC simulation with the best fit Ma are shown in Fig. 
with all five energy regions combined. The best fit values 
for the free parameters in the fit are summarized in Tab. 

In the fit, there is a strong correlation between the nor- 
malization, nonQE/QE, and the two-track to one-track 
migration, but different ways of constraining the param- 
eters do not affect the fit value of Ma very much. The 
migration is expressed such that 0.90 means 10% of the 
two-track events in each E rec and Q 2 ec bin should be 
moved to the corresponding one-track bin; this param- 
eter is being pulled beyond its 5% gaussian constraint 
for both samples. For example, fixing the 2-track to 1- 
track migration at 0.95, the one-sigma edge of the 5% 
error, yields normalization = 0.975, nonQE/QE = 1.26, 
and proton rescattering 1.20. However, Ma is unchanged 
while the chisquare rises to 261 for 234 degrees of free- 
dom. 

Finally, one further comment about the nonQE/QE 
parameter. The higher fit value for Af^ E =1.20 GeV 
causes an increase of about 10% in the absolute QE 



cross section, relative to the default value of mi =1.10 
GeV. If a good fit requires maintaining a similar relative 
non-QE cross-section, then a corresponding increase in 
nonQE/QE parameter and decrease in the absolute nor- 
malization is required. This appears as a portion of the 
nonQE/QE = 1.30 fit value and normalization. 



parameter 


fit value error 


M A (GeV) 


1.20 


0.09 


$(0.5 to 1.0 GeV) 


1.02 


0.25 


$(1.0 to 1.5 GeV) 


1.00 


fixed 


$(1.5 to 2.0 GeV) 


0.80 


0.09 


$(2.0 to 2.5 GeV) 


0.93 


0.08 


$(> 2.5 GeV) 


1.09 


0.11 


Normalization 


0.96 


0.09 


nonQE/QE 


1.30 


0.17 


K2K-I 2tk -> ltk 


0.94 


0.02 


K2K-IIa 2tk -> ltk 


0.94 


0.03 


Proton Rescattering 


1.14 


0.18 



TABLE IV: Best fit values for the parameters in the fit. The 
errors given are from the fit only. The error for Ma rises to 
±0.12 GeV when the other systematic effects are included 



B. Consistency checks 

A test for consistency is to vary the low Q 2 cut and 
compare the results. In Fig. [SJ the best fit Ma is shown 
with different minimum Q 2 ec . The error bars include an 
estimate of statistical errors only; however, the data for 
each point in this figure are correlated with the other 
points. The extra error bar shows the total error. When 
no cut is applied (and no coherent pion), the fit value is 
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Ma = 1-27 ± 0.12 GeV, where the statistical error is less, 
but a large additional systematic error of ± 0.07 GeV is 
assigned (but not shown in Fig. |HJ) due to uncertainty 
in the amount of Pauli blocking. When the systematic 
errors are considered, these results are consistent with 
the fit value. 
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FIG. 8: Fit values obtained for different values of the low Q 2 
cut. Only statistical errors are shown. The horizontal line 
is the combined best fit. The vertical line is the systematic 
errors. 



A second check is to consider the fit values for the 
Q 2 distribution at each energy, shown in Fig. [§J This 
uses the best fit values for the flux for all energies except 
the one being tested while the chisquare, and therefore 
the shape fit, is computed only for the energy bins in 
question. This is necessary because of the significant mi- 
gration from true energy (where the flux parameter is 
applied) to reconstructed energy used in the fit, espe- 
cially for the non-QE background. There are different 
systematic effects, and this result should not be consid- 
ered a measurement, but rather a consistency test. How- 
ever results for each energy are also consistent with the 
combined result. 

Finally, we fit for the value of Ma for the K2K-I and 
K2K-IIa data sets separately, and obtain the values 1.12 
± 0.12 (x 2 150/127 dof) and 1.25 ± 0.18 GeV (x 2 = 
109/101 dof) respectively. The primary difference be- 
tween the two data sets are the presence or absence of 
the lead glass detector and the accompanying uncertainty 
in the muon momentum scale. Also, the muon thresh- 
olds are somewhat different, and the acceptance model 
for short, second tracks in SciFi is calibrated separately 
for the two data sets, which could show up in the 2-track 
to 1-track migration. 



C. Systematic uncertainties 

Here these systematic errors are discussed in detail. 
The largest contributions to the systematic error, sum- 
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FIG. 9: Fit values obtained separately from the shape of the 
Q 2 distribution for each neutrino energy. The horizontal line 
indicates the combined best fit value. 



marized in Tab.lVl are the uncertainty in the muon mo- 
mentum scale, and the normalization and uncertainty in 
the flux for each energy region. Other, smaller contribu- 
tions include the shape of the non-QE background, the 
non-QE/QE ratio, and the two-track to one-track mi- 
gration. A final, interesting source of uncertainty comes 
from nuclear effects, though it contributes only a small 
amount to this analysis. The statistical error is esti- 
mated by setting all the other parameters in the fit to 
their best fit values and determining the resulting error 
in Ma, though there is a further statistical effect in the 
normalization parameter. 



Sources of uncertainty 


Error in MA 


Muon momentum scale 


0.07 


Relative Flux and normalization 


0.06 


M A 1-7T 


0.03 


nonQE/QE 


0.03 


Proton rescattering 


0.03 


Statistics 


0.03 


Total 


0.12 



TABLE V: The calculation of the total error. Errors smaller 
than 0.03 are not listed. The total value takes into account 
the correlations among those errors that are parameters in 
the Ma fit; the others are added to that total in quadrature. 



1. Muon momentum scale 

The muon momentum appears directly, and indirectly 
via E v , in the calculation of the value of Q 2 ec for each 
event. The uncertain absolute scale for this momentum, 
as modeled in the detector Monte Carlo simulation, will 
cause the MC prediction for the shape of the Q 2 ec distri- 
bution to be more or less compressed and affect the Ma 
measurement. As an example, a ± 1% error in the mo- 
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mentum scale gives a =F 0.05 GeV error in the fit value 
for Ma- Approximately =p 0.01 GeV of this error can 
be attributed to shifting a small number of events up or 
down one E r v ec bin. The other =p 0.04 GeV is from the 
calculation of the reconstructed Q 2 itself. The central 
value for the muon momentum scale is determined from 
the neutrino energy spectrum analysis, described in Sec. 
Ill and reference [17], while the error from that analysis 
is propagated to the Ma analysis as described below. 

Because the muon momentum is measured using its 
range in the detector, the uncertainty for the overall mo- 
mentum may come from any of the pieces of the detec- 
tor: SciFi, LG, SciBar, or the MRD. In this analysis, 
we model this uncertainty by assigning it to just two of 
these. The first is the uncertainty in the density of the 
lead-glass detector and therefore the energy loss experi- 
enced by the muon passing through it. The second piece 
is a scaling factor for part of the muon momentum cal- 
culated from the range in the MRD detector. For both 
pieces, we determine the central value of the momentum 
shift and the error from the neutrino data. 

The density of the lead glass, which is incorporated 
into the geometry description in our Monte Carlo simu- 
lation, is determined from a beam test and is uncertain 
by 5%. We have modeled the effect of this uncertainty 
and made a reweighting table that modifies the MC p M 
and 6^ distribution. This uncertainty could give rise to 
a 2% error in the total momentum for a typical K2K-I 
event. In the neutrino energy spectrum measurement, 
this is a parameter in the fit and good agreement with 
the data is found with a value that is 0.98 ± 0.013 times 
the density obtained from the beam test; the neutrino 
data provide the stronger constraint. This central value 
is used in the Ma analysis. 

Likewise, we measure a shift in the muon momentum 
scale for the Muon Range Detector of 0.976 ± 0.007 using 
the energy spectrum measurement procedure. When the 
K2K-I and K2K-IIa data are fit separately, we obtain a 
consistent result for this parameter, despite the presence 
of the lead glass detector in the former. This is assigned 
as an error for the MRD portion of the muon range, but 
it actually arises from a combination of factors includ- 
ing the material assay for the MRD and SciFi (about 
1%), the simulation of muon energy loss in GEANT [24| 
(about 1%) and the uncertainty in the intrinsic muon mo- 
mentum from the neutrino interaction MC (about 0.5%). 
Again, we find the neutrino data produce a good central 
value and a tighter constraint than taking the individual 
errors in quadrature. Though these errors actually come 
from all portions of the muon track, we find no signifi- 
cant difference in the analysis if this factor is obtained 
from and applied to the whole track momentum, instead 
of the MRD portion only. 

Because the Ma fit and the energy spectrum analy- 
sis use the same neutrino data, it is possible that the 
uncertain value for Ma itself is affecting the fit values 
for the MRD muon momentum scale when that value is 
obtained from the spectrum measurement. Our default 



Monte Carlo assumes Ma = 1.1 GeV. An uncertainty 
in this value of ± 0.20 GeV corresponds to an error of ± 
0.01 in the fit value of the momentum scale. This is taken 
as an additional uncertainty when this parameter is used 
to determine Ma- Also, there is a correlation between 
the lead-glass density error and the MRD momentum er- 
ror. When all of these effects are combined, the resulting 
error in M A is ±0.07 GeV. 



2. Flux for each energy region and normalization 

A significant uncertainty in Ma arises because the rel- 
ative neutrino flux for each energy region and the overall 
normalization parameters are unconstrained parameters 
in the fit. In this way we are fitting the shape of the Q 2 
distribution in each energy region separately, regardless 
of the errors in the incident neutrino flux. The contribu- 
tion to the total error is estimated by fixing the other free 
parameters in the fit and reading the resulting error in 
Ma, which is from these parameters and statistics only. 

The overall normalization contributes more to the error 
than the uncertain relative normalization. This is esti- 
mated by further constraining the relative flux so that 
only the normalization and Ma are free. The overall 
normalization is correlated with Ma because Ma affects 
both the relative size of and the shape of the QE cross- 
section. Different combinations of Ma and normalization 
will give a reasonable chisquare when compared with the 
data, and the error due to this parameter, more than the 
others, would be reduced with increased data statistics, 
even with no further constraints. This last result is con- 
firmed using MC samples of various sizes as if they were 
data, to study the effect of statistics of the data sample. 

We do have a constraint on the relative flux for each 
energy region from the neutrino oscillation measurement 
This measurement is done using data from all the 
near detectors, not just SciFi. This information is not 
completely independent of this analysis because it shares 
some of the same data set, but a different analysis tech- 
nique, and several other data sets from the other near 
detectors. We get a consistent result Ma = 1.13 ± 0.12 
GeV when this constraint is used. 



3. nonQE/QE parameter and two-track to one-track 
migration 

The nonQE/QE scaling ratio is also a free parame- 
ter in the axial- mass fit. There is no constraint on this 
parameter for this analysis, thou gh o ther estimates find 
that it is uncertain by 5 to 10% 17]. The two-track to 
one-track migration parameter is highly correlated with 
nonQE/QE, and when these two are combined, they con- 
tribute a total error of 0.03 GeV to the fit value of Ma- 
As before, this is obtained by fixing all the other param- 
eters such that the resulting error in Ma is the combina- 
tion of these and the statistical error only. 
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4- Non quasi-elastic background shape 

Single pion events from the production and decay of 
the A and other resonances in the nucleus are the largest 
background to the QE samples in this analysis. These 
events are described by a calculation that includes a sim- 
ilar axial mass parameter which affects the shape of the 
Q 2 distribution. If the value used to model the single pion 
background is different, that will affect the fit value ob- 
tained for the quasi-elastic events. Our calculation takes 
Mjf = 1.1 ± 0.1 GeV. This contributes an uncertainty of 
± 0.03 to result for M® E , and is estimated by generating 
a second complete MC sample with M A V =1.2 GeV. 

Other contributions to the nonQE background are 
deep inelastic scattering and coherent pion production. 
For the former, we have evaluated the uncertainty by re- 
moving the Bodek-Yang correction, and find no effect. 
We also consider the case where charged-current coher- 
ent pion events are produced which has only +0.01 GeV 
effect for the Q 2 > 0.2 cut used in the standard analysis, 
but increases the fit value by 0.10 GeV when we fit the 
entire Q 2 range. 

5. Nuclear effects 

Other interesting sources of uncertainty are the effects 
of the nucleus on the cross section and the Q 2 distribu- 
tion, primarily from the nucleon momentum distribution. 
The effects arc small relative to the other uncertainties 
described above because the minimum Q 2 ec cut elimi- 
nates the data where these errors are most significant. 
These effects will be of interest for future precision ex- 
periments and as models of neutrino-nucleus interactions 
become more sophisticated. We present a description of 
these effects for the uniform Fermi gas model, in this case 
from the calculation in [2^, |2(j . The three effects are de- 
scribed below and summarized in Fig. ^|for a 1.0 GeV 
neutrino. It is the ratio in this figure that is the basis for 
R(E U ,Q 2 ) in Eq. El 

The main uncertainty is the amount of Pauli blocking 
that should be applied both to the quasi-elastic and also 
the single pion background. Within the context of the 
Fermi gas model, this can be estimated by assuming a 
different kf. 215 and 235 in addition to the default value 
of 225 GeV/c. The effects of this uncertainty on Ma 
do not appear with the Q 2 > 0.2 requirement used in 
this analysis, but are as much as 5% for fits that use the 
lowest Q 2 events. 

At the upper end of the Q 2 distribution, the quasi- 
elastic cross section has a kinematic cut off whose loca- 
tion depends on the incident neutrino energy. The mo- 
mentum distribution in a nucleus smears this step, giving 
a tail to the distribution. These high Q 2 interactions pro- 
duce muons that do not reach the MRD because they are 
at high angle or their momentum is too low, so this has 
no effect on the present analysis. 

The momentum distribution will modify the shape of 



the Q 2 distribution through the middle region between 
the two effects described in the preceding paragraphs. 
The slope of the middle region in the second plot in Fig. 
1101 is approximately 0.017 (GeV/c)~ 2 . There is also an 
overall suppression of the cross section of 2%. The uncer- 
tainty represented by the change in slope can be propa- 
gated to the Ma analysis by modifying R(E,Q 2 ) in the fit. 
The resulting uncertainty in Ma is ±0.01 GeV, negligible 
compared to the other uncertainties in this analysis. 

A final uncertainty from the nuclear model is the nu- 
cleon interaction energy. For our Fermi gas model, this 
takes the form of an effective binding energy -27 ± 3 
MeV, and is the energy given up to the recoil proton from 
the nucleus. This affects the outgoing muon momentum 
and would contribute ± 0.02 GeV error to Ma, but this 
is naturally included by the free muon momentum scale 
parameter in this analysis. 

These uncertainties are also used to estimate the ef- 
fect of the 21.8% aluminum that makes up the fiducial 
mass. The neutrino-aluminum interactions are taken to 
have the same cross section per nucleon and the same 
kinematics as for oxygen. A higher kf appropriate for 
aluminum only has an effect in the Pauli blocked region. 
The increased effective binding energy is equivalent to a 
shift in p^ of about 3 MeV for this fraction of the inter- 
actions, and thus is negligible for the whole sample. 



D. Effect of the new vector form factors 

The basic method used to measure the axial vector 
mass here is the same as for previous measurements, 
listed in Tab. EH but since that time there have been 
improved measurements for the shape of the vector form 
factors from electron scattering experiments. Changing 
the shape of the contribution of the vector form factors 
affects the fit shape of the axial-vector form factor. We 
continue to assume Fa has a dipole form; future high- 
precision neutrino experiments may be sensitive to subtle 
differences from a Fa dipole shape. 

Our results assume the updated parameterizations 
of Bosted @- We have evaluated one other recent 
parametrization and found the Ma result differs by 
only 0.01 GeV. This is true even considering the discrep- 
ancy between the polarization transfer measurement and 
the Rosenbluth separation measurement, described in Q 
with further references. For our analysis, we use the pa- 
rameterization of given by [33j . 

In order to allow comparison with the previous results, 
we have repeated the analysis with the same modified 
dipole approximations used by [JIEE1 who follow Olsson 
et al. 34]. We find that these old parameterizations 
produce a value that is 1.23, roughly 0.03 higher. When 
only a pure dipole is used, the fit value is 1.27 GeV. 
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FIG. 10: Effect of the nucleon momentum on the shape of the Q distribution. The comparison is between the free nucleon 
and a uniform Fermi gas model. The effect of Pauli blocking is seen at low Q 2 , the tail of the momentum distribution at high 
Q 2 , an overall suppression, and a slight change in the slope in the middle region. The calculated quasi-elastic cross sections for 
1.0 GeV neutrinos on oxygen are on the left, and the ratio (Fermi gas)/ (free neutron). 



Experiment Pub. Date 


Target 


Method 


M A 


Error 


comment 


ANL [6J 


1982 


D 


12' Bubble Chamber 


1.00 


± 0.05 




FNAL [5] 


1983 


D 


15' Bubble Chamber 


1.05 


+0.12 - 0.16 




BNL [4] 


1990 


D 


7' Bubble Chamber 


1.07 +0.040 -0.045 




CERN [2J] 


1977 


CF 3 Br 


GGM Bubble Chamber 


0.94 


± 0.17 




CERN [28] 


1979 


CF 3 Br, C 3 H 8 


GGM Bubble Chamber 


0.94 


± 0.05 




SKAT [22] 


1990 


CF 3 Br 


Bubble Chamber 


1.05 


± 0.14 


(") 


SKAT [29] 


1990 


CF 3 Br 


Bubble Chamber 


0.79 


± 0.20 


(*) 


BNL [30] 


1969 


Fe 


Segmented Tracker 


1.05 


± 0.20 




BNL [21] 


1987 


HC, Al 


Segmented Tracker 


1.06 


± 0.05 


elastic scattering 


BNL [32] 


1988 


HC, Al 


Segmented Tracker 


1.09 


± 0.04 


(?) 


K2K SciFi 


this expt. 


EhO, Al 


Segmented Tracker 


1.27 


± 0.12 


dipole form factors 



TABLE VI: Results from other experiments, grouped first by target nucleus, then by publication date. Where separate values 
are given for Ma extracted from the shape of da/dQ 2 only, that is the value included in this table. All the data are for the 
neutrino quasi-elastic reaction (y n — ► fj,~ p) except for two which also took data with anti-neutrino (p p — *• fi + n), one of which 
studied neutral current (elastic) scattering, noted in the table. For better comparison with other experiments, the K2K SciFi 
result is the one analyzed with dipole vector form factors. 



E. Comparison with other experiments 



tions for a selection of these experiments. 



This is the first measurement of the axial vector mass 
using neutrino interactions with oxygen targets, but 
there have been many previous measurements with a va- 
riety of other target nuclei. The experiments in Tab. IVII 
have hundreds or thousands of events from neutrino or 
anti-neutrinos with energies of a few GeV. The system- 
atic errors they report are dominated by uncertainties 
in the neutrino flux, calculation of nuclear effects, and 
subtraction of non-quasielastic backgrounds. 



One problem with comparing the results in Tab. IVII 
is that the older analyses used not only different as- 
sumptions about the vector form factors, but also dif- 
ferent backgrounds and other physical constants such as 
F A (q 2 =0). The results given here are the published val- 
ues; however, the authors of [3j have made some effort 
to reproduce and then update all of the analysis assump- 



VI. CONCLUSION 

We have made the first measurement of axial vector 
form factor using neutrino interactions on an oxygen tar- 
get. We find an axial vector mass Ma — 1-20 ± 0.12 GeV 
gives the best agreement with the data, if we assume a 
dipole form for Fa- This analysis includes the updated 
(non-dipole) vector form factors obtained from electron 
scattering experiments. In order to better compare with 
previous experiments, an alternate result using only pure 
dipole vector form factors is M A = 1.27 ± 0.12 GeV. We 
note that this analysis is very sensitive to the absolute 
muon momentum scale. We have also studied the details 
of the nucleon momentum distribution for oxygen on this 
analysis and find only a small effect on the shape of the 
Q 2 distribution for Q 2 > 0.2 (GeV/c) 2 . 
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